Parametric Resonance in the Drift Motion of an Ionic Bubble in near critical Ar Gas 
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The drift mobility of the OJT ion in dense Argon gas near the liquid-vapor critical point has been 
measured as a function of the density. Near T c the zero-field density-normalized ion mobility fioN 
shows a deep minimum at a density smaller than N c . This phenomenon was previously intepreted 
as the result of the formation of a correlated cluster of Argon atoms around the ion because of the 
strong electrostriction exerted by the ion on the highly polarizable and compressible gas. We now 
suggest that a possible alternative explanation is related to the onset of a parametric resonance of 
a bubble surrounding the ion. At resonance a large amount of energy is dissipated by sound waves 
in addition to viscous dissipation processes, resulting in the large mobility drop observed. 
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The transport properties of negative ions in dense rare 
gases and liquids have been recently subject of renewed 
interest because they can give information on how the 
microscopic structure of the fluid around the ion is mod- 
ified by the ion-atom interaction Jl],^. Moreover, the 
possibility to greatly change with relative ease the gas 
density N close to the critical point of the liquid-vapor 
transition gives the experimenters the unique opportu- 
nity to investigate the transition from the kinetic to the 
hydrodynamic transport regime Q| . 

The drift mobility of in dense Ne gas Q or in liquid 
Xc Q can be satisfactorily explained in terms of hydro- 
dynamics if ions are assumed to be surrounded by a com- 
plex structure. The competition between short-range re- 
pulsive exchange forces between the weakly bound outer 
electron of the ion and the electrons of the atoms and 
the long-range polarization attraction of the ion with the 
atoms leads to the formation of a microcavity around the 
ion, which is in turn surrounded by a strong density en- 
hancement [|J||. 

A self-consistent field model has been developed to 
compute the structure shape p|. The size of the com- 
plex is taken as effective hydrodynamic radius R of the 
ion and its drift mobility is calculated by means of the 
Stokes formula fiQ = e/Q-KrjR, where r\ is the viscosity. 
The agreement with mobility data in liquid Xe 0] is 
quite good, but the data in Ne gas near T c ||] are almost 
quantitatively reproduced only for N > N c . In particu- 
lar, the model does not reproduce even qualitatively the 
little dimple in the density-normalized mobility /j,qN for 
N « 0.7N C . In this Letter we show that this general fea- 
ture of the mobility || can be explained by assuming that 
under certain termodynamic conditions a further dissi- 
pation mechanism becomes active in addition to viscous 
processes, namely sound wave emission by oscillations of 
the structure surrounding the ion. 

Since the influence of electrostriction on the structure 
formation depends on the gas polarizability, we have car- 
ried out mobility measurements in Ar gas because 



its polarizability is much larger than that of Ne. We used 
the pulsed photoinjection technique as for Ne [Q] and He 
[3|. Details of the experiment are reported in literature 
[3,9. A small bunch of electrons is injected into the gas 
by irradiating a photocathode by means of a short UV 
light pulse. Electrons are captured by O2 impurities (in a 
concentration of 10— 100 p. p.m.) forming stable ions 
p0| . They drift under the action of an external electric 
field E towards the anode inducing a detectable current. 
The analysis of its time dependence allows the determi- 
nation of the drift time and, hence, of the drift velocity 
v d p3|. Finally, the mobility is calculated as fi = vd/E. 
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FIG. 1. Experimental zero-field density-normalized mobil- 
ity fi N vs. A7iV c for T = 151.5 K (T/T w 1.005). 

In Fig|l] we show the measured zero-field density- 
normalized mobility (IqN vs. the reduced density N/N c 
for T= 151.5 K(T/T C » 1.005) (JV C = 8.08 atoms • nm" 3 
T c = 151. 7K). hqN shows a very deep minimum for 
N/N c ps 0.76. The Stokes formula with constant radius 
reproduces the data, at most, for N/N c > 1.3 @. The 
data are heuristically explained by assuming that a clus- 
ter of Ar atoms forms around the ion as a consequence of 
electrostriction M . This large structure interacts hydro- 
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dynamically with the gas. The effective hydro-dynamic 
radius depends on the local properties of the fluid. The 
good agreement with the data for N/N c > 0.5 is obtained 
by introducing into the effective radius an adjustable con- 
tribution proportional to the size of the strongly cor- 
related cluster. Thus, the hydrodynamic radius shows 
a complicated ad hoc density dependence, whose mean- 
ing is not easy to grasp. The reason might be that this 
model focuses only on the structure determined by the 
long-range part of the ion-atom interaction potential and 
does not take into account the presence of the microcav- 
ity closely surrounding the ion. Therefore, it neglects 
the possibility that the cavity may oscillate. In this case 
the energy of the moving complex, ion plus structure or, 
briefly, ionic bubble, can be dissipated by sound wave ra- 
diation in addition to viscous processes. 

The present problem bears many affinities with that of 
the dynamics of electron bubble formation in liquid He 
fjJIl, the adiabatic process where the repulsive electron- 
atom interaction displaces a large number of atoms far 
away from the electron which is localized within the fluid 
dilation, or the electron solvation dynamics in water, 
where cavity contraction is induced by long-range attrac- 
tive polarization interactions fjl|. These problems have 
been treated in terms of a hydrodynamic model of cavity 
expansion jl2| or collapse |l3| . The dynamics of the cav- 
ity boundary is determined by the flow of the surround- 
ing solvent. This is assumed to be spherically symmetric 
with velocity in the radial direction. In absence of energy 
dissipation and by assuming the incompressibility of the 
solvent, the cavity boundary velocity U is described by 
the Rayleigh-Plesset (RP) equation 
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where R is the time dependent cavity radius, p is the 
gas mass density, p e is the external pressure, and p(R) is 
the pressure on the cavity boundary, expressed in terms 
of the free energy change F in the process: p(R) = 
— (dF/dR) /4ttR 2 . The equilibrium radius is obtained 
by minimizing F, but in absence of dissipation the sys- 
tem will not equilibrate. The cavity oscillates around 
the equilibrium value of the radius. Energy dissipation 
causes the damping of the oscillations. Within the given 
approximations, the first passage time, i.e., the time re- 
quired for the cavity radius to expand (contract) from its 
initial value to its maximum (minimum) is taken as the 
lower bound for the relaxation timescale of the process. 
Viscosity and emission of sound waves drive the system 
towards the equilibrium Jl2|,[l3| ■ 

It is obvious that also the drift motion of the ion is 
affected if energy is dissipated also by emission of sound 
waves. We therefore address the issue of the motion of the 
ionic bubble in dense Ar gas, considering the oscillations 
of the bubble boundary. We assume a spherical cavity 



with well-defined boundary of radius R(t) with equilib- 
rium value Rq. The main contribution to the restoring 
force is due to the surface tension. Although a gas has no 
surface tension, the enhanced density region around the 
cavity acts as a well-defined interface. Oscillations are 
initiated by collisions with energetic Ar atoms. Owing to 
the bubble smallness and to the asimmetry induced by 
the average motion driven by the external electric field, 
these collisions are not very frequent and occur in a spa- 
tially non homogeneous way. The equation of motion of 
the bubble boundary is then given by the RP equation, 
and, beside the impulsive nature of the driving force, is 
the same equation used for the bubble dynamics in the 
sonoluminescence experiments [|l5| . 

So, for small distortion of the spherical shape a solution 
for the bubble radius is sought in the form Rq + a n Y n , 
where Y n is a spherical harmonic of degree n, and a n 
are the distortion amplitude coefficients. As bubble os- 
cillations are initiated by collisions with energetic Ar 
atoms, this phenomenon can be basically described by 
a kicked rotator model, with the following equation of 
motion for the simple normalized periodic case: a" + 
raj, — Kf(a n ) X^m=o — mr ) = 0' wnere m is an inte- 
ger, r the damping constant and r the period between 
two kicks. Imposing a stochastic behaviour to the po- 
tential function, we can model collision events occurring 
with a gaussian distribution, e.g. in the quantum kicked 
rotator and electron localization problem |m| . All prop- 
erties of this equation are well known results of nonlinear 
dynamical system theory, which assure, under certain as- 
sumptions, the existence of a universal route to chaos. In 
the limit of small forcing the dynamics of the distortion 
amplitude can be cast in the form of the Mathieu-Hill 
(MH) like equation p) 
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with b n oc R^ 2 a n {t). uj 2 n = {loq/uj) 2 is the square of the 
ratio between the natural frequency of the bubble, luq, 
and the excitation frequency u>. The natural frequency is 
given by ui 2 = f3 n a/pRl, where /3„ = (n-l)(n+l)(n+2), 
a is the surface tension, and p is the mass density of the 
gas jl5[ . Primes denote differentiation with respect to the 
dimensionless time t ~ ujt. Assuming that the bubble is 
empty, the surface tension can be calculated by using the 
parachoric formula [|l6| a — (P7V) 4 , where N is the gas 
number density and P is a constant. For Argon P w 
1.39 x 1CT 29 J 1 / 4 m 5 / 2 . The term f m = 2n(n + 2)ri/puR 2 l 
is a damping coefficient related to the viscosity rj whose 
values are found in literature p7y . Of the forcing term 
only the first Fourier component of amplitude e m has 
been retained. 

For a dissipation-free system the Mathieu Hill equa- 
tion is known to give origin to parametric instability, 
when deviations from the spherical shape accumulate 
over many oscillation cycles. Floquet's theorem states 



2 



that solutions of Eqj^ take the form jl8| 
b n (t) = e^P n (t) 



(3) 



where [if is the Floquet's index and P„ is periodic. w m 
and e m span a plane geometrically divided into stability 
and instability regions |l{J. We focus on the n — 2 mode. 
In the instability regions fi 2 > and the envelope of the 
solutions to the MH equation grows exponentially. In 
the physical systems at hand this amplitude cannot grow 
indefinitely because dissipation stabilizes the surface dy- 
namics. In the stability regions where Re /if — and 
ImjjLf = 2u) m the solutions are periodic with angular fre- 
quency ujq with a small modulation at angular frequency 
2luq. In this case the system experiences parametric res- 
onance. In the limit of small forcing (e m — > 0), in ab- 
sence of dissipation (£ m — * 0) and close to the resonance 
{uj m ~ 1), the non-zero solution for the Floquet's index 
is fif w 2iu> m — 2iu)$/u) Jl9| | - We now show that Im/if 
is related to the measured mobility /^o. From the expres- 
sion of the bubble eigenfrequency, using the parachoric 
formula, we can write 



(4) 



where g = (N A f3 n P 4 /MR 3 Q /2 ) 1 / 2 is a constant. N A is 
the Avogadro's number and M is the atomic weight 
of Argon, ui is the excitation frequency. Therefore, 
under the hypotesis of collisional induced oscillations, 
uj = lo co u/1 = 2nv co u/l, where v co ii is the collision fre- 
quency and / defines the order of a suitable subharmonic. 
In the Knudsen regime the mobility of a heavy ion of ra- 
dius Rq scattered off light particles of mass m is 
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By inserting this result into Eq.|| we get 



Im/i/ cx VN(noN) 
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where fi^N is the experimentally measured zero-field 
density-normalized mobility. Therefore, yN (j,qN can be 
used to make a stability analysis of the bubble boundary 
motion. 

In the real system, the ideal conditions of small forc- 
ing and absence of dissipation are not met and a para- 
metric resonance is characterized by a minimum of the 
effective Floquet index, approximated by Re [if, which 
is shown to be related to y/~N fioN. The minimum ob- 
served for N/N c sa 0.76 (N sa 6.2 atoms • nm 3 ) is due 
to parametric resonance of the ionic bubble. We follow 



a heuristic line of reasoning. Once excited, the bubble 
starts oscillating with characteristic frequency Im{/i/} 
and growing amplitude given by A(t) = exp (Re{[if}t). 
The growth of the amplitude is limited because collision 
processes can absorb the oscillation energy. Therefore, 
if T is the time interval during which the bubble oscil- 
lates, Re{/i/}T ~ O(l). On the other hand, according to 
calculations in the cavitation model [^3|, the oscillations 
last for a few cycles, hence Im{/i/}T ~ 0(1). Then, we 
have Re{^/} ~ Im{^/} and, owing to Eq]^, y/N(n N) 
is also related to the real part of the Floquet's index. 
The Ansatz must be then verified by a direct numerical 
stability analysis of Eq.^. 

The minimum of /iqN occurs at a density N/N c « 0.76 
where the resonance condition u — ujq is met. In Fig.^ 



we plot vs. iV both ujq oc N 3 / 2 and lo 
approximately. 
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FIG. 2. Excitation frequency ui and eigenfrequency uo of 
the oscillating ionic bubble. 

The two curves intersect for the right density value if 
I w 2 11 . So, the excitation frequency must be w 2 11 times 
lower than the average collision frequency given by Eqj^. 
The factor 2 11 can be explained assuming that only colli- 
sions with Ar atoms in the tail of the Maxwell-Boltzmann 
energy distribution function, namely those with average 
kinetic energy in excess of ik^T, are energetic enough to 
initiate bubble oscillations. These energetic collisions are 
less frequent than average by the right factor. 

We have numerically integrated the MH equation 
(Eq.||) by assuming an equilibrium radius Rq 10 A, 
and have carried out the usual analysis of the stability 
bands. Usually, the Floquet's index is plotted vs. uo^. 
However, to compare the experimental data with the re- 
sults of numerical analysis we have converted uj m to N by 
means of Eqns. ||and||. In Fig. || we plot (/.i N)y/N and 
Refif as a function of the gas density N. The Floquet's 
index shows a number of small bands, but also a very 
deep minimum at the same N of the experimental /iqN. 
This minimum is therefore associated with a strong para- 
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metric resonance of the ionic bubble. The larger width 
of the experimental data can be explained by considering 
the fact that the experimental result is an average over 
a distribution of ionic bubble radii. At the condition 
of parametric resonance the amplitude of oscillation can 
reach substantial values and the kinetic energy gained by 
the external electric field can be efficiently dissipated by 
emission of sound waves. 



of the bubble boundary. We finally note that in this 
problem there is a moving interface between media of 
different polarizability, crossed by the strong electrical 
field of the ion. Therefore, an experiment could be de- 
signed to detect the expected quantum radiation emitted 
as a dynamic Casimir effect, that some authors consider 
the physical cause of sonoluminescence |23]. 
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FIG. 3. Floquet's index (if obtained by numerical integrat- 
ing the MH equation with Ro = 10 A and t m = 1100 com- 
pared to the experimental data (iqN\/~N vs. N/N c . 

Close to the critical point the sound velocity c is mini- 
mum (w 190m/s [p0|) and sound dissipation is favoured. 
The sound intensity irradiated by an oscillating bubble 
in the long wavelength limit is [£lj / = 2ttpcll> 2 Rq. At the 
density of the minimum fioN the intensity irradiated per 
cycle is in the tens of meV range and a relevant part of 
the energy gained by the ion from the electric field is dis- 
sipated by means of this process in addition to the usual 
viscous processes. At higher N the stiffness of the bubble 
surface becomes too large and the oscillations cannot be 
initiated so easily as in the resonance region. Hence, at 
high N the contribution of the sound to energy dissipa- 
tion is limited and (Iq is determined by the viscosity. This 
is why the high— N region /iqN can be well described by 
the Stokes formula. On the contrary, at much lower N 
bubbles cannot form and fJ®N is determined by dissipa- 
tion processes other than sound emission. 

In the region where stable bubbles exist a further dis- 
sipation mechanism, though probably very small, should 
be considered. Within the bubble the ion undergoes a 
chaotic motion bouncing back and forth from the inner 
bubble wall. Molecular Dynamics studies p3] have shown 
a vibration of the ion within the cavity. The ion should 
therefore behave as an emitting antenna of characteristic 
frequency uo e w 2irv th /2R , where v t h = (3k B T/mi) 1/2 
is the ion thermal velocity and mi its mass. At the ex- 
periment temperature uj e w 10 rads. Moreover, this 
radiation should be modulated by the slower oscillation 
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